\(\int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx\) [509]

   Optimal result
   Rubi [A] (verified)
   Mathematica [C] (warning: unable to verify)
   Maple [B] (verified)
   Fricas [A] (verification not implemented)
   Sympy [F(-1)]
   Maxima [C] (verification not implemented)
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 45, antiderivative size = 237 \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\frac {\sqrt {2} (A-B+C) \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right )}{\sqrt {a} d}+\frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \]

[Out]

(A-B+C)*arctan(1/2*sin(d*x+c)*a^(1/2)*2^(1/2)/cos(d*x+c)^(1/2)/(a+a*cos(d*x+c))^(1/2))*2^(1/2)/d/a^(1/2)+2/7*A
*sin(d*x+c)/d/cos(d*x+c)^(7/2)/(a+a*cos(d*x+c))^(1/2)-2/35*(A-7*B)*sin(d*x+c)/d/cos(d*x+c)^(5/2)/(a+a*cos(d*x+
c))^(1/2)+2/105*(31*A-7*B+35*C)*sin(d*x+c)/d/cos(d*x+c)^(3/2)/(a+a*cos(d*x+c))^(1/2)-2/105*(43*A-91*B+35*C)*si
n(d*x+c)/d/cos(d*x+c)^(1/2)/(a+a*cos(d*x+c))^(1/2)

Rubi [A] (verified)

Time = 0.87 (sec) , antiderivative size = 237, normalized size of antiderivative = 1.00, number of steps used = 7, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.111, Rules used = {3122, 3063, 12, 2861, 211} \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\frac {\sqrt {2} (A-B+C) \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a \cos (c+d x)+a}}\right )}{\sqrt {a} d}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a \cos (c+d x)+a}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a \cos (c+d x)+a}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a \cos (c+d x)+a}}+\frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a \cos (c+d x)+a}} \]

[In]

Int[(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)/(Cos[c + d*x]^(9/2)*Sqrt[a + a*Cos[c + d*x]]),x]

[Out]

(Sqrt[2]*(A - B + C)*ArcTan[(Sqrt[a]*Sin[c + d*x])/(Sqrt[2]*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Cos[c + d*x]])])/(Sq
rt[a]*d) + (2*A*Sin[c + d*x])/(7*d*Cos[c + d*x]^(7/2)*Sqrt[a + a*Cos[c + d*x]]) - (2*(A - 7*B)*Sin[c + d*x])/(
35*d*Cos[c + d*x]^(5/2)*Sqrt[a + a*Cos[c + d*x]]) + (2*(31*A - 7*B + 35*C)*Sin[c + d*x])/(105*d*Cos[c + d*x]^(
3/2)*Sqrt[a + a*Cos[c + d*x]]) - (2*(43*A - 91*B + 35*C)*Sin[c + d*x])/(105*d*Sqrt[Cos[c + d*x]]*Sqrt[a + a*Co
s[c + d*x]])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 2861

Int[1/(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]]), x_Symbol] :> D
ist[-2*(a/f), Subst[Int[1/(2*b^2 - (a*c - b*d)*x^2), x], x, b*(Cos[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c +
 d*Sin[e + f*x]]))], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 -
 d^2, 0]

Rule 3063

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])*((c_.) + (d_.)*sin[(e_
.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[(B*c - A*d)*Cos[e + f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x]
)^(n + 1)/(f*(n + 1)*(c^2 - d^2))), x] + Dist[1/(b*(n + 1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^m*(c + d*Sin
[e + f*x])^(n + 1)*Simp[A*(a*d*m + b*c*(n + 1)) - B*(a*c*m + b*d*(n + 1)) + b*(B*c - A*d)*(m + n + 2)*Sin[e +
f*x], x], x], x] /; FreeQ[{a, b, c, d, e, f, A, B, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2
 - d^2, 0] && LtQ[n, -1] && (IntegerQ[n] || EqQ[m + 1/2, 0])

Rule 3122

Int[((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_)*((A_.) + (B_.)*s
in[(e_.) + (f_.)*(x_)] + (C_.)*sin[(e_.) + (f_.)*(x_)]^2), x_Symbol] :> Simp[(-(c^2*C - B*c*d + A*d^2))*Cos[e
+ f*x]*(a + b*Sin[e + f*x])^m*((c + d*Sin[e + f*x])^(n + 1)/(d*f*(n + 1)*(c^2 - d^2))), x] + Dist[1/(b*d*(n +
1)*(c^2 - d^2)), Int[(a + b*Sin[e + f*x])^m*(c + d*Sin[e + f*x])^(n + 1)*Simp[A*d*(a*d*m + b*c*(n + 1)) + (c*C
 - B*d)*(a*c*m + b*d*(n + 1)) + b*(d*(B*c - A*d)*(m + n + 2) - C*(c^2*(m + 1) + d^2*(n + 1)))*Sin[e + f*x], x]
, x], x] /; FreeQ[{a, b, c, d, e, f, A, B, C, m}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^
2, 0] &&  !LtQ[m, -2^(-1)] && (LtQ[n, -1] || EqQ[m + n + 2, 0])

Rubi steps \begin{align*} \text {integral}& = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 \int \frac {-\frac {1}{2} a (A-7 B)+\frac {1}{2} a (6 A+7 C) \cos (c+d x)}{\cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{7 a} \\ & = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {4 \int \frac {\frac {1}{4} a^2 (31 A-7 B+35 C)-a^2 (A-7 B) \cos (c+d x)}{\cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{35 a^2} \\ & = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {8 \int \frac {-\frac {1}{8} a^3 (43 A-91 B+35 C)+\frac {1}{4} a^3 (31 A-7 B+35 C) \cos (c+d x)}{\cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx}{105 a^3} \\ & = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}+\frac {16 \int \frac {105 a^4 (A-B+C)}{16 \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \, dx}{105 a^4} \\ & = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}+(A-B+C) \int \frac {1}{\sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \, dx \\ & = \frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}-\frac {(2 a (A-B+C)) \text {Subst}\left (\int \frac {1}{2 a^2+a x^2} \, dx,x,-\frac {a \sin (c+d x)}{\sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right )}{d} \\ & = \frac {\sqrt {2} (A-B+C) \arctan \left (\frac {\sqrt {a} \sin (c+d x)}{\sqrt {2} \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}}\right )}{\sqrt {a} d}+\frac {2 A \sin (c+d x)}{7 d \cos ^{\frac {7}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (A-7 B) \sin (c+d x)}{35 d \cos ^{\frac {5}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}+\frac {2 (31 A-7 B+35 C) \sin (c+d x)}{105 d \cos ^{\frac {3}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}}-\frac {2 (43 A-91 B+35 C) \sin (c+d x)}{105 d \sqrt {\cos (c+d x)} \sqrt {a+a \cos (c+d x)}} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains higher order function than in optimal. Order 5 vs. order 3 in optimal.

Time = 9.46 (sec) , antiderivative size = 2797, normalized size of antiderivative = 11.80 \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\text {Result too large to show} \]

[In]

Integrate[(A + B*Cos[c + d*x] + C*Cos[c + d*x]^2)/(Cos[c + d*x]^(9/2)*Sqrt[a + a*Cos[c + d*x]]),x]

[Out]

(4*B*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(7*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(7/2)
) - (4*C*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(7*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(
7/2)) + (24*B*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(35*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/2 + (d*x)/2
]^2)^(5/2)) + (4*C*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(35*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/2 + (d
*x)/2]^2)^(5/2)) + (32*B*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(35*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/
2 + (d*x)/2]^2)^(3/2)) + (16*C*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(105*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2
*Sin[c/2 + (d*x)/2]^2)^(3/2)) + (64*B*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(35*d*Sqrt[a*(1 + Cos[c + d*x])]*
Sqrt[1 - 2*Sin[c/2 + (d*x)/2]^2]) + (32*C*Cos[c/2 + (d*x)/2]*Sin[c/2 + (d*x)/2])/(105*d*Sqrt[a*(1 + Cos[c + d*
x])]*Sqrt[1 - 2*Sin[c/2 + (d*x)/2]^2]) + (2*(A - B + C)*Cot[c/2 + (d*x)/2]*Csc[c/2 + (d*x)/2]^8*(363825*Sin[c/
2 + (d*x)/2]^2 - 4729725*Sin[c/2 + (d*x)/2]^4 + 26785605*Sin[c/2 + (d*x)/2]^6 - 86790165*Sin[c/2 + (d*x)/2]^8
+ 177677808*Sin[c/2 + (d*x)/2]^10 - 239283044*Sin[c/2 + (d*x)/2]^12 + 52080*Hypergeometric2F1[2, 11/2, 13/2, S
in[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12 + 560*HypergeometricPFQ[{2, 2, 2, 2,
11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12 + 213120160*
Sin[c/2 + (d*x)/2]^14 - 168280*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2
]^2)]*Sin[c/2 + (d*x)/2]^14 - 2240*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2
/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^14 - 121497024*Sin[c/2 + (d*x)/2]^16 + 212520*Hypergeometri
c2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^16 + 3360*Hypergeom
etricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d
*x)/2]^16 + 40125184*Sin[c/2 + (d*x)/2]^18 - 124320*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1
+ 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^18 - 2240*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2},
 Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^18 - 5840384*Sin[c/2 + (d*x)/2]^20 + 2
8000*Hypergeometric2F1[2, 11/2, 13/2, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^2
0 + 560*HypergeometricPFQ[{2, 2, 2, 2, 11/2}, {1, 1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]
^2)]*Sin[c/2 + (d*x)/2]^20 + 363825*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sqrt[Sin
[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 5336100*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 +
 (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^2*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 34636140*ArcTan
h[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^4*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1
 + 2*Sin[c/2 + (d*x)/2]^2)] - 131060160*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[
c/2 + (d*x)/2]^6*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 320535600*ArcTanh[Sqrt[Sin[c/2 + (
d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^8*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*
x)/2]^2)] - 530671680*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^10*
Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 604296000*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2
*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^12*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 4689
48480*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^14*Sqrt[Sin[c/2 + (
d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] + 237726720*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)
/2]^2)]]*Sin[c/2 + (d*x)/2]^16*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 70963200*ArcTanh[Sqr
t[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 + (d*x)/2]^18*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2
*Sin[c/2 + (d*x)/2]^2)] + 9461760*ArcTanh[Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]]*Sin[c/2 +
(d*x)/2]^20*Sqrt[Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)] - 1120*Cos[(c + d*x)/2]^6*HypergeometricP
FQ[{2, 2, 2, 11/2}, {1, 1, 13/2}, Sin[c/2 + (d*x)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12*(-
6 + 5*Sin[c/2 + (d*x)/2]^2) + 280*Cos[(c + d*x)/2]^4*HypergeometricPFQ[{2, 2, 11/2}, {1, 13/2}, Sin[c/2 + (d*x
)/2]^2/(-1 + 2*Sin[c/2 + (d*x)/2]^2)]*Sin[c/2 + (d*x)/2]^12*(103 - 164*Sin[c/2 + (d*x)/2]^2 + 70*Sin[c/2 + (d*
x)/2]^4)))/(40425*d*Sqrt[a*(1 + Cos[c + d*x])]*(1 - 2*Sin[c/2 + (d*x)/2]^2)^(9/2)*(-1 + 2*Sin[c/2 + (d*x)/2]^2
))

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(474\) vs. \(2(202)=404\).

Time = 13.17 (sec) , antiderivative size = 475, normalized size of antiderivative = 2.00

method result size
default \(-\frac {\left (105 A \sqrt {2}\, \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )-105 B \sqrt {2}\, \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )+105 \sqrt {2}\, \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, C \left (\cos ^{4}\left (d x +c \right )\right ) \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+105 A \left (\cos ^{3}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \sqrt {2}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )-105 B \left (\cos ^{3}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \sqrt {2}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+105 C \left (\cos ^{3}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \sqrt {2}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+86 A \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )-182 B \sin \left (d x +c \right ) \left (\cos ^{3}\left (d x +c \right )\right )+70 C \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )-62 A \sin \left (d x +c \right ) \left (\cos ^{2}\left (d x +c \right )\right )+14 B \sin \left (d x +c \right ) \left (\cos ^{2}\left (d x +c \right )\right )-70 C \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right )+6 A \sin \left (d x +c \right ) \cos \left (d x +c \right )-42 B \sin \left (d x +c \right ) \cos \left (d x +c \right )-30 A \sin \left (d x +c \right )\right ) \sqrt {a \left (1+\cos \left (d x +c \right )\right )}}{105 d a \left (1+\cos \left (d x +c \right )\right ) \cos \left (d x +c \right )^{\frac {7}{2}}}\) \(475\)
parts \(-\frac {A \left (105 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right ) \left (\cos ^{4}\left (d x +c \right )\right )+43 \sqrt {2}\, \left (\cos ^{3}\left (d x +c \right )\right ) \sin \left (d x +c \right )+105 \left (\cos ^{3}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )-31 \sqrt {2}\, \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right )+3 \sin \left (d x +c \right ) \cos \left (d x +c \right ) \sqrt {2}-15 \sqrt {2}\, \sin \left (d x +c \right )\right ) \sqrt {a \left (1+\cos \left (d x +c \right )\right )}\, \sqrt {2}}{105 d \left (1+\cos \left (d x +c \right )\right ) \cos \left (d x +c \right )^{\frac {7}{2}} a}+\frac {B \left (15 \left (\cos ^{3}\left (d x +c \right )\right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+13 \sqrt {2}\, \left (\cos ^{2}\left (d x +c \right )\right ) \sin \left (d x +c \right )+15 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \left (\cos ^{2}\left (d x +c \right )\right ) \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )-\sin \left (d x +c \right ) \cos \left (d x +c \right ) \sqrt {2}+3 \sqrt {2}\, \sin \left (d x +c \right )\right ) \sqrt {a \left (1+\cos \left (d x +c \right )\right )}\, \sqrt {2}}{15 d \left (1+\cos \left (d x +c \right )\right ) \cos \left (d x +c \right )^{\frac {5}{2}} a}-\frac {C \left (3 \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \left (\cos ^{2}\left (d x +c \right )\right ) \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )+\sin \left (d x +c \right ) \cos \left (d x +c \right ) \sqrt {2}+3 \cos \left (d x +c \right ) \sqrt {\frac {\cos \left (d x +c \right )}{1+\cos \left (d x +c \right )}}\, \arcsin \left (\cot \left (d x +c \right )-\csc \left (d x +c \right )\right )-\sqrt {2}\, \sin \left (d x +c \right )\right ) \sqrt {a \left (1+\cos \left (d x +c \right )\right )}\, \sqrt {2}}{3 d \left (1+\cos \left (d x +c \right )\right ) \cos \left (d x +c \right )^{\frac {3}{2}} a}\) \(539\)

[In]

int((A+B*cos(d*x+c)+C*cos(d*x+c)^2)/cos(d*x+c)^(9/2)/(a+cos(d*x+c)*a)^(1/2),x,method=_RETURNVERBOSE)

[Out]

-1/105/d*(105*A*2^(1/2)*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*arcsin(cot(d*x+c)-csc(d*x+c))*cos(d*x+c)^4-105*B*2^(
1/2)*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*arcsin(cot(d*x+c)-csc(d*x+c))*cos(d*x+c)^4+105*2^(1/2)*(cos(d*x+c)/(1+c
os(d*x+c)))^(1/2)*C*cos(d*x+c)^4*arcsin(cot(d*x+c)-csc(d*x+c))+105*A*cos(d*x+c)^3*(cos(d*x+c)/(1+cos(d*x+c)))^
(1/2)*2^(1/2)*arcsin(cot(d*x+c)-csc(d*x+c))-105*B*cos(d*x+c)^3*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*2^(1/2)*arcsi
n(cot(d*x+c)-csc(d*x+c))+105*C*cos(d*x+c)^3*(cos(d*x+c)/(1+cos(d*x+c)))^(1/2)*2^(1/2)*arcsin(cot(d*x+c)-csc(d*
x+c))+86*A*cos(d*x+c)^3*sin(d*x+c)-182*B*sin(d*x+c)*cos(d*x+c)^3+70*C*cos(d*x+c)^3*sin(d*x+c)-62*A*sin(d*x+c)*
cos(d*x+c)^2+14*B*sin(d*x+c)*cos(d*x+c)^2-70*C*cos(d*x+c)^2*sin(d*x+c)+6*A*sin(d*x+c)*cos(d*x+c)-42*B*sin(d*x+
c)*cos(d*x+c)-30*A*sin(d*x+c))*(a*(1+cos(d*x+c)))^(1/2)/a/(1+cos(d*x+c))/cos(d*x+c)^(7/2)

Fricas [A] (verification not implemented)

none

Time = 0.31 (sec) , antiderivative size = 205, normalized size of antiderivative = 0.86 \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=-\frac {2 \, {\left ({\left (43 \, A - 91 \, B + 35 \, C\right )} \cos \left (d x + c\right )^{3} - {\left (31 \, A - 7 \, B + 35 \, C\right )} \cos \left (d x + c\right )^{2} + 3 \, {\left (A - 7 \, B\right )} \cos \left (d x + c\right ) - 15 \, A\right )} \sqrt {a \cos \left (d x + c\right ) + a} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right ) - \frac {105 \, \sqrt {2} {\left ({\left (A - B + C\right )} a \cos \left (d x + c\right )^{5} + {\left (A - B + C\right )} a \cos \left (d x + c\right )^{4}\right )} \arctan \left (\frac {\sqrt {2} \sqrt {a \cos \left (d x + c\right ) + a} \sqrt {\cos \left (d x + c\right )} \sin \left (d x + c\right )}{2 \, {\left (\cos \left (d x + c\right )^{2} + \cos \left (d x + c\right )\right )} \sqrt {a}}\right )}{\sqrt {a}}}{105 \, {\left (a d \cos \left (d x + c\right )^{5} + a d \cos \left (d x + c\right )^{4}\right )}} \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)/cos(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="fricas")

[Out]

-1/105*(2*((43*A - 91*B + 35*C)*cos(d*x + c)^3 - (31*A - 7*B + 35*C)*cos(d*x + c)^2 + 3*(A - 7*B)*cos(d*x + c)
 - 15*A)*sqrt(a*cos(d*x + c) + a)*sqrt(cos(d*x + c))*sin(d*x + c) - 105*sqrt(2)*((A - B + C)*a*cos(d*x + c)^5
+ (A - B + C)*a*cos(d*x + c)^4)*arctan(1/2*sqrt(2)*sqrt(a*cos(d*x + c) + a)*sqrt(cos(d*x + c))*sin(d*x + c)/((
cos(d*x + c)^2 + cos(d*x + c))*sqrt(a)))/sqrt(a))/(a*d*cos(d*x + c)^5 + a*d*cos(d*x + c)^4)

Sympy [F(-1)]

Timed out. \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\text {Timed out} \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)**2)/cos(d*x+c)**(9/2)/(a+a*cos(d*x+c))**(1/2),x)

[Out]

Timed out

Maxima [C] (verification not implemented)

Result contains complex when optimal does not.

Time = 1.39 (sec) , antiderivative size = 3203, normalized size of antiderivative = 13.51 \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\text {Too large to display} \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)/cos(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="maxima")

[Out]

1/105*((105*(sqrt(2)*cos(2*d*x + 2*c)^4 + sqrt(2)*sin(2*d*x + 2*c)^4 + 4*sqrt(2)*cos(2*d*x + 2*c)^3 + 2*(sqrt(
2)*cos(2*d*x + 2*c)^2 + 2*sqrt(2)*cos(2*d*x + 2*c) + sqrt(2))*sin(2*d*x + 2*c)^2 + 6*sqrt(2)*cos(2*d*x + 2*c)^
2 + 4*sqrt(2)*cos(2*d*x + 2*c) + sqrt(2))*sqrt(a)*arctan2(((abs(e^(I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 + sin(
d*x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d*x +
 c)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) + 1)*sin(d*x + c)^2 + 6*cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1/4)*
sin(1/2*arctan2(2*(cos(d*x + c) - 1)*sin(d*x + c)/abs(e^(I*d*x + I*c) + 1)^2, (abs(e^(I*d*x + I*c) + 1)^2 + co
s(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)/abs(e^(I*d*x + I*c) + 1)^2)) + sin(d*x + c))/abs(e^(I*d*x
+ I*c) + 1), ((abs(e^(I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 + sin(d*x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)
^2 - 2*cos(d*x + c) + 1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d*x + c)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) +
1)*sin(d*x + c)^2 + 6*cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1/4)*sqrt(a)*cos(1/2*arctan2(2*(cos(d*x + c) - 1)*
sin(d*x + c)/abs(e^(I*d*x + I*c) + 1)^2, (abs(e^(I*d*x + I*c) + 1)^2 + cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos
(d*x + c) + 1)/abs(e^(I*d*x + I*c) + 1)^2)) + sqrt(a)*cos(d*x + c) - sqrt(a))/(sqrt(a)*abs(e^(I*d*x + I*c) + 1
))) - 2*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(3/4)*(43*(cos(2*d*x + 2*c)^2*sin(d
*x + c) + sin(2*d*x + 2*c)^2*sin(d*x + c) + 2*cos(2*d*x + 2*c)*sin(d*x + c) + sin(d*x + c))*cos(1/2*arctan2(si
n(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) + 108*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))*sin(d*x
+ c) - 12*(9*cos(d*x + c) - 7)*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - ((43*cos(d*x + c) -
105)*cos(2*d*x + 2*c)^2 + (43*cos(d*x + c) - 105)*sin(2*d*x + 2*c)^2 + 2*(43*cos(d*x + c) - 105)*cos(2*d*x + 2
*c) + 43*cos(d*x + c) - 105)*sin(1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))*sqrt(a) - 4*(cos(2*d*x
+ 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*(12*(14*sin(2*d*x + 2*c) - 5*sin(d*x + c))*cos(7
/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) + 37*(cos(2*d*x + 2*c)^2*sin(d*x + c) + sin(2*d*x + 2*c)^2
*sin(d*x + c) + 2*cos(2*d*x + 2*c)*sin(d*x + c) + sin(d*x + c))*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x +
2*c) + 1)) - 12*(14*cos(2*d*x + 2*c) - 5*cos(d*x + c) + 9)*sin(7/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c)
+ 1)) - ((37*cos(d*x + c) - 35)*cos(2*d*x + 2*c)^2 + (37*cos(d*x + c) - 35)*sin(2*d*x + 2*c)^2 + 2*(37*cos(d*x
 + c) - 35)*cos(2*d*x + 2*c) + 37*cos(d*x + c) - 35)*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)))
*sqrt(a))*A/(a*cos(2*d*x + 2*c)^4 + a*sin(2*d*x + 2*c)^4 + 4*a*cos(2*d*x + 2*c)^3 + 6*a*cos(2*d*x + 2*c)^2 + 2
*(a*cos(2*d*x + 2*c)^2 + 2*a*cos(2*d*x + 2*c) + a)*sin(2*d*x + 2*c)^2 + 4*a*cos(2*d*x + 2*c) + a) + 35*(3*(sqr
t(2)*cos(2*d*x + 2*c)^2 + sqrt(2)*sin(2*d*x + 2*c)^2 + 2*sqrt(2)*cos(2*d*x + 2*c) + sqrt(2))*arctan2(((abs(e^(
I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 + sin(d*x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) +
1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d*x + c)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) + 1)*sin(d*x + c)^2 + 6*
cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1/4)*sin(1/2*arctan2(2*(cos(d*x + c) - 1)*sin(d*x + c)/abs(e^(I*d*x + I*
c) + 1)^2, (abs(e^(I*d*x + I*c) + 1)^2 + cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)/abs(e^(I*d*x +
I*c) + 1)^2)) + sin(d*x + c))/abs(e^(I*d*x + I*c) + 1), ((abs(e^(I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 + sin(d*
x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d*x + c
)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) + 1)*sin(d*x + c)^2 + 6*cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1/4)*sq
rt(a)*cos(1/2*arctan2(2*(cos(d*x + c) - 1)*sin(d*x + c)/abs(e^(I*d*x + I*c) + 1)^2, (abs(e^(I*d*x + I*c) + 1)^
2 + cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)/abs(e^(I*d*x + I*c) + 1)^2)) + sqrt(a)*cos(d*x + c)
- sqrt(a))/(sqrt(a)*abs(e^(I*d*x + I*c) + 1))) - 2*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*
c) + 1)^(3/4)*(cos(1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))*sin(d*x + c) - (cos(d*x + c) - 3)*sin(
1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))) - 4*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d
*x + 2*c) + 1)^(1/4)*(cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))*sin(d*x + c) - (cos(d*x + c) +
1)*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))))*C/((cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*
cos(2*d*x + 2*c) + 1)*sqrt(a)) - 7*(15*(sqrt(2)*cos(2*d*x + 2*c)^2 + sqrt(2)*sin(2*d*x + 2*c)^2 + 2*sqrt(2)*co
s(2*d*x + 2*c) + sqrt(2))*(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos(2*d*x + 2*c) + 1)^(1/4)*arctan2(((a
bs(e^(I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 + sin(d*x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x +
 c) + 1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d*x + c)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) + 1)*sin(d*x + c)^
2 + 6*cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1/4)*sin(1/2*arctan2(2*(cos(d*x + c) - 1)*sin(d*x + c)/abs(e^(I*d*
x + I*c) + 1)^2, (abs(e^(I*d*x + I*c) + 1)^2 + cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)/abs(e^(I*
d*x + I*c) + 1)^2)) + sin(d*x + c))/abs(e^(I*d*x + I*c) + 1), ((abs(e^(I*d*x + I*c) + 1)^4 + cos(d*x + c)^4 +
sin(d*x + c)^4 + 2*(cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)*abs(e^(I*d*x + I*c) + 1)^2 - 4*cos(d
*x + c)^3 + 2*(cos(d*x + c)^2 - 2*cos(d*x + c) + 1)*sin(d*x + c)^2 + 6*cos(d*x + c)^2 - 4*cos(d*x + c) + 1)^(1
/4)*sqrt(a)*cos(1/2*arctan2(2*(cos(d*x + c) - 1)*sin(d*x + c)/abs(e^(I*d*x + I*c) + 1)^2, (abs(e^(I*d*x + I*c)
 + 1)^2 + cos(d*x + c)^2 - sin(d*x + c)^2 - 2*cos(d*x + c) + 1)/abs(e^(I*d*x + I*c) + 1)^2)) + sqrt(a)*cos(d*x
 + c) - sqrt(a))/(sqrt(a)*abs(e^(I*d*x + I*c) + 1))) - 26*(cos(2*d*x + 2*c)^2*sin(d*x + c) + sin(2*d*x + 2*c)^
2*sin(d*x + c) + 2*cos(2*d*x + 2*c)*sin(d*x + c) + sin(d*x + c))*cos(1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x +
 2*c) + 1)) + 24*cos(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))*sin(d*x + c) - 24*(cos(d*x + c) - 1)
*sin(5/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) + 2*((13*cos(d*x + c) - 15)*cos(2*d*x + 2*c)^2 + (13
*cos(d*x + c) - 15)*sin(2*d*x + 2*c)^2 + 2*(13*cos(d*x + c) - 15)*cos(2*d*x + 2*c) + 13*cos(d*x + c) - 15)*sin
(1/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1)) - 4*sqrt(cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 + 2*cos
(2*d*x + 2*c) + 1)*(7*cos(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))*sin(d*x + c) - (7*cos(d*x + c)
- 5)*sin(3/2*arctan2(sin(2*d*x + 2*c), cos(2*d*x + 2*c) + 1))))*B/((cos(2*d*x + 2*c)^2 + sin(2*d*x + 2*c)^2 +
2*cos(2*d*x + 2*c) + 1)^(5/4)*sqrt(a)))/d

Giac [F]

\[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\int { \frac {C \cos \left (d x + c\right )^{2} + B \cos \left (d x + c\right ) + A}{\sqrt {a \cos \left (d x + c\right ) + a} \cos \left (d x + c\right )^{\frac {9}{2}}} \,d x } \]

[In]

integrate((A+B*cos(d*x+c)+C*cos(d*x+c)^2)/cos(d*x+c)^(9/2)/(a+a*cos(d*x+c))^(1/2),x, algorithm="giac")

[Out]

integrate((C*cos(d*x + c)^2 + B*cos(d*x + c) + A)/(sqrt(a*cos(d*x + c) + a)*cos(d*x + c)^(9/2)), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {A+B \cos (c+d x)+C \cos ^2(c+d x)}{\cos ^{\frac {9}{2}}(c+d x) \sqrt {a+a \cos (c+d x)}} \, dx=\int \frac {C\,{\cos \left (c+d\,x\right )}^2+B\,\cos \left (c+d\,x\right )+A}{{\cos \left (c+d\,x\right )}^{9/2}\,\sqrt {a+a\,\cos \left (c+d\,x\right )}} \,d x \]

[In]

int((A + B*cos(c + d*x) + C*cos(c + d*x)^2)/(cos(c + d*x)^(9/2)*(a + a*cos(c + d*x))^(1/2)),x)

[Out]

int((A + B*cos(c + d*x) + C*cos(c + d*x)^2)/(cos(c + d*x)^(9/2)*(a + a*cos(c + d*x))^(1/2)), x)